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Abstract 

We consider the gluon fusion production cross section of a scalar Higgs boson at NLO QCD in the 
SM and in the MSSM. We implement the calculation in the POWHEG approach, and match the NLO- 
QCD results with the PYTHIA and HERWIG QCD parton showers. We discuss a few representative 
scenarios in the SM and MSSM parameter spaces, with emphasis on the fermion and squark mass 
effects on the Higgs boson distributions. 
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1 Introduction 



Understanding the mechanism that leads to the breaking of the electroweak symmetry and that is 
responsible for the generation of the mass of the elementary particles is one of the major challenges 
of high energy physics. The search for the Higgs boson(s) is currently under way at the Tevatron and 
at the LHC, and limits on the Higgs mass spectrum have already been set [HE]- This search requires 
an accurate control of all the Higgs production and decay mechanisms, including the effects due to 
radiative corrections [3]. 

In the Standard Model (SM) the gluon fusion process [4] is the dominant Higgs production mecha- 
nism both at the Tevatron and at the LHC. The total cross section receives very large next-to- leading 
order (NLO) QCD corrections, which were first computed in ref. [5] in the so-called heavy-quark ef- 
fective theory (HQET), i.e. including only the top-quark contributions in the limit writ — > oo. Later 
calculations [3 El QU] retained the exact dependence on the masses of the top and bottom quarks 
running in the loops. The next-to-next-to-leading order (NNLO) QCD corrections are also large, and 
have been computed in the HQET in ref. [11]. The finite-top-mass effects at NNLO QCD have been 
studied in ref. [12] and found to be small. The resummation to all orders of soft gluon radiation 
has been studied in refs. \13\ 114] , Leading third-order (NNNLO) QCD terms have been discussed 
in ref. |15j . The role of electroweak (EW) corrections has been discussed in refs. |16t \T7\ \18\ 119] . 
The impact of mixed QCD-EW corrections has been discussed in ref. [20]. The residual uncertainty 
on the total cross section depends mainly on the uncomputed higher-order QCD effects and on the 
uncertainties that affect the parton distribution functions (PDF) of the proton [2"T j l2"2" j 13]. 

The Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) consists of two SU (2) 
doublets, Hi and H2, whose relative contribution to electroweak symmetry breaking is determined by 
the ratio of vacuum expectation values of their neutral components, tan f3 = V2/V1. The spectrum of 
physical Higgs bosons is richer than in the SM, consisting of two neutral CP-even bosons, h and H, one 
neutral CP-odd boson, A, and two charged bosons, H ± . The couplings of the MSSM Higgs bosons 
to matter fermions differ from those of the SM Higgs, and they can be considerably enhanced (or 
suppressed) depending on tan/3. As in the SM, gluon fusion is one of the most important production 
mechanisms for the neutral Higgs bosons, whose couplings to the gluons are mediated by top and 
bottom quarks and their supersymmetric partners, the stop and sbottom squarks. 

In the MSSM, the cross section for Higgs boson production in gluon fusion is currently known at 
the NLO. The contributions arising from diagrams with quarks and gluons can be obtained from the 
corresponding SM results [Q El O [TU] with an appropriate rescaling of the Higgs-quark couplings. 
The contributions arising from diagrams with squarks and gluons were first computed under the 
approximation of vanishing Higgs mass in ref. [23] , and the full Higgs-mass dependence was included 
in later calculations [8j [91 QUI [21] • The contributions of two-loop diagrams involving top, stop and 
gluino to both scalar and pseudoscalar Higgs production were computed under the approximation 
of vanishing Higgs mass in refs. [25\ 126] . whose results were later confirmed and cast in a compact 
analytic form in refs. [27\ 128] . 
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The approximation of vanishing Higgs mass can provide reasonably accurate results as long as 
the Higgs mass is well below the threshold for creation of the massive particles running in the loops. 
For the production of the lightest scalar Higgs, this condition does apply to the two-loop diagrams 
involving top, stop and gluino, but it obviously does not apply to the corresponding diagrams involving 
the bottom quark, whose contribution can be relevant for large values of tan f3. In turn, the masses 
of the heaviest scalar and of the pseudoscalar might very well approach (or exceed) the threshold for 
creation of top quarks or even of squarks. Unfortunately, retaining the full dependence on the Higgs 
mass in the quark-squark-gluino contributions has proved a rather daunting task. A calculation based 
on a combination of analytic and numerical methods was presented in ref. |29j (see also ref. |30j). 
but neither explicit analytic results nor a public computer code have been made available so far. 
However, ref. |31| presented an approximate evaluation of the bottom-sbottom-gluino contributions 
to scalar production, based on an asymptotic expansion in the large super symmetric masses that is 
valid up to and including terms of 0(m%/rnV), 0(mt,/MsusY) and 0{m 2 z /Mg USY ), where denotes 
a Higgs boson mass and Msusy denotes a generic superparticle mass. An independent calculation of 
the bottom-sbottom-gluino contributions, restricted to the limit of a degenerate superparticle mass 
spectrum, was also presented in ref. |32| . confirming the results of ref. [31J. More recently, ref. [28J 
presented an evaluation of the quark-squark-gluino contributions to pseudoscalar production that is 
also based on an asymptotic expansion in the large super symmetric masses, but does not assume any 
hierarchy between the pseudoscalar mass and the quark mass, thus covering both the top-stop-gluino 
and bottom-sbottom gluino cases. 

The total cross section, without acceptance cuts, provides important information about the Higgs 
boson production rate. On the other hand, especially at the LHC, it is very likely that a Higgs boson is 
produced in association with a jet, generating a transverse momentum p!f of the Higgs boson. The first 
studies on Higgs+jet final states in the SM were performed in ref. [33], using the real-parton emission 
amplitudes that enter the calculation of the NLO-QCD corrections to the inclusive Higgs production 
cross section. The NLO corrections to the Higgs+jet final state at large transverse momentum of the 
Higgs boson were subsequently studied in ref. [34] . and the resummation of all the logarithmically 
enhanced terms, matched with the NLO calculation of the Higgs+jet final state, was discussed in 
ref. [35J. The impact on the Higgs+jet final state of the NLO-EW corrections and of the finite masses 
of the particles running in the loops was discussed in ref. [36] • In the MSSM, the production of a 
neutral Higgs boson in association with one jet was discussed in refs. [57| 158]. 

A different set of observables are the differential distributions of the Higgs boson, inclusive over 
QCD radiation. For the SM case, results were presented in ref. [39J at NLO-QCD and in ref. [30] at 
NNLO-QCD. The Higgs boson transverse momentum spectrum, including the NLO-QCD corrections 
matched with the resummation of next-to-next-to-leading logarithmic (NNLL) enhanced terms, has 
been studied in refs. [3TJ[321[33]. In the MSSM a study of the Higgs distributions, at NLO-QCD, was 
discussed in ref. [32J. 

If a new scalar particle is discovered at the Tevatron or at the LHC, a major question will be to 
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determine whether it is a Higgs boson and, in that case, whether it belongs to the particle spectrum 
of the SM, of the MSSM or of any other model. An example could be represented by a MSSM Higgs 
boson whose production cross section is close to the production cross section for a SM Higgs boson of 
equal mass. In this case, an accurate study of the differential distributions involving the Higgs boson 
might shed some light on the underlying model. 

A precise analysis of the experimental data requires the use of NLO-QCD results merged with 
the description of initial-state multiple gluon emission via a QCD parton shower (PS) Monte Carlo 
(MC) such as HERWIG [2] or PYTHIA [35]. However, the merging of NLO-QCD matrix elements with 
PS faces the problem of avoiding double counting, as addressed in refs. [46j|47]. The POWHEG method 
[48j allows to systematically merge NLO calculations with vetoed PS, avoiding double counting and 
preserving the NLO accuracy of the calculation. The procedure can be implemented using a set of 
tools and results available in the so-called POWHEG BOX The latter provides a general framework 
that exploits the universal nature of initial-state collinear divergences and the factorization property 
of soft radiation to automatize the subtraction of all the soft and/or collinear divergent terms from 
the NLO matrix elements of an arbitrary process. The POWHEG method does not rely on the details 
of the shower MC and, by construction, guarantees an accuracy at NLO + leading logarithmic (LL) 
QCD. In ref. [51] it has been shown that, with an appropriate choice of scale for the strong coupling 
constant, the merging procedure can also reproduce the next-to-leading logarithmic (NLL) terms. 

At present, no code exists that merges the NLO-QCD results for the gluon fusion process with a 
QCD PS, retaining the exact dependence on the Higgs mass and on the masses of the particles running 
in the loops. Two implementations of the NLO-QCD results merged with a PS are availabl^]: the one 
in MCONLO |50] and the one in the POWHEG BOX framework [51]. However, both implementations are 
limited to the SM case, and beyond LO they only include results computed in the HQET. Conversely, 
the codes HIGLU [53J and iHixs [M] contain the full dependence on the Higgs and quark masses up to 
NLO, and HIGLU also allows to include the full squark-gluon contributions from ref. [24J, but neither 
code is matched to a shower MC. Recently, a step toward the inclusion of the finite-quark- mass effects 
in PS was taken in ref. [55J, where parton- level events for Higgs production accompanied by zero, 
one or two partons are generated with matrix elements computed in the HQET, and then, before 
being passed to the PS, they are re- weighted by the ratio of the exact one- loop amplitudes over the 
approximate ones. This procedure is equivalent to generating events directly with the exact one-loop 
amplitudes, yet it is much faster. 

We aim to provide a code that fills the remaining gap, using matrix elements that include the 
dependence on the masses, both in the SM and in the MSSM, properly matched to an external shower 
MC. For the SM case, we use NLO matrix elements with full dependence on the Higgs, top and bottom 
masses. For the MSSM case, we use matrix elements with exact dependence on the quark, squark and 
Higgs masses in the contributions of real-parton emission. For the two-loop virtual contributions, the 
approximation of vanishing Higgs mass is employed in the diagrams involving superpartners, while 

1 These implementations are also available as a subprocess of HERWIG++ |52] , 
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the rest is computed exactly. 

The plan of the paper is as follows: in section [2] we describe the basic features of the POWHEG 
implementation of gg — > (j>; in section [3] we discuss our SM implementation with exact dependence on 
the fermion masses, presenting a numerical analysis valid for an on-shell Higgs; section [5] is devoted 
to analyzing the MSSM case; finally, in section [5] we draw our conclusions. 



2 POWHEG implementation of <?<?—>■ (j) 

In this section we briefly discuss the implementation of the gluon-fusion Higgs production process in 
the POWHEG BOX framework, following closely ref. [51] (see also ref. [56]). We fix the notation keeping 
the discussion at a general level, without referring to a specific model. In the next sections the formulae 
presented below will be specialized to the SM and MSSM cases. 

The generation of the hardest emission is done in POWHEG according to the following formula: 

da = d§a | A + A ($ 1)PT ) d$rad j 

+ R m (*1> $ rad) d$i<f$ rad . (1) 

q 

In the equation above the variables $i = (M 2 , Y) denote the invariant mass squared and the rapidity 
of the Higgs boson, which describe the kinematics of the Born (i.e., lowest-order) process gg — > (j). 
The variables $ ra d = (£ 5 2/j </>) describe the kinematics of the additional final-state parton in the real 
emission processes. In particular, denoting by k' 2 the momentum of the final-state parton in the 
partonic center-of-mass frame, or 

k' 2 = k' 2 ° (1, sin 9sm<j), sin 9 cos (f>, cos 9), (2) 

we have 

4° = ^, y = cos9, (3) 

where s is the partonic center-of-mass energy squared. 

The factor B(<&i) in eq. ([T]) is related to the total cross section computed at NLO in QCD. 
It contains the value of the differential cross section, for a given configuration of the Born final 
state variables, integrated over the radiation variables. The integral of this quantity on d&i without 
acceptance cuts yields the total cross section. This factor is responsible for the correct NLO-QCD 
normalization of the result, and is computed in the initialization phase using the real and virtual 
NLO-QCD corrections. 

The terms within curly brackets in eq. ([1]) describe the real emission spectrum of an extra parton: 
the first term is the probability of not emitting any parton with transverse momentum larger than 
a cutoff py m , while the second term is the probability of not emitting any parton with transverse 
momentum larger than a given value px times the probability of emitting a parton with transverse 
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momentum equal to px- The sum of the two terms fully describes the probability of having either zero 
or one additional parton in the final state. The probability of non-emission of a parton with transverse 
momentum kx larger than is obtained using the POWHEG Sudakov form factor 

#(*l,$rad) 



A($i,p T ) = exp 



d$rad- 



-9{hr ~Pt) 



(4) 



Finally, the last term in eq. ([!]) describes the effect of the qq — > 4>g channel, which has been kept 
apart in the generation of the first hard emission because it does not factorize into the Born cross 
section times an emission factor. 

We now discuss the various terms appearing in eq. (H|) in more detail. We have: 



= B gg (^) + V gg (^) 

+ J d^rad | Rgg (*1, $rad) + ^ [^79 (*!' $rad ) + ^<?9 (*!' $rad ) 



+ c. r. , 



where 



V w (*i)4rc> 
7^ ffff (#i,$ rad )£ ffff , 

^ OT (*l»*«d)^OTJ 

^ ?s (#i,$ rad )£ ?s , 



(5) 



(6) 
(7) 
(8) 
(9) 
(10) 



with £ a b the luminosity for the partons a and b. In eq. © "c. r." denotes the collinear remnants 
multiplied by the relevant parton luminosity. The remnants are the finite leftovers after the subtraction 
of the initial-state collinear singularities into the parton distribution function is performed, and their 
explicit expressions are given in eqs. (2.36), (2.37) and (3.7)-(3.10) of ref. [51] . 

The function B gg (<&i) in eq. © represents the squared matrix element of the Born contribution 
to the process, averaged over colors and helicities of the incoming gluons, and multiplied by the flux 
factor 1/(2M 2 ). It is given by 



H 



(11) 



256\/2 7r 2 

where % is the form factor for the coupling of the Higgs boson with two gluons, whose explicit form 
depends on the particle content of the model considered and will be detailed in the following sections. 
It is decomposed in one- and two loop parts as 



U 



, a s aj2e , i ,,2 , 



TT 



+ 0(ai) . 



The regularized two-loop virtual contributions are contained in 

r2 / ,,2 \ 



V, 



99 



($0 = 2s 

7T 



(12) 



(13) 
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In the equation above, fiR and fip are the renormalization and factorization scale, respectively, Ca = 
N c (N c being the number of colors), and /3q = (11 Ca — 2 Nf) /6 (Nf being the number of active flavors) 
is the one-loop beta function of the strong coupling. 

The hatted functions IZij in eqs. ([8))- (fT0]) are the Frixione, Kunst and Signer (57J EE] infrared- 
subtracted counterparts of IZij 

^(*i,*«d) = j \l Q) [(iz^) + (rb)J } [(1 " 2 /2 )^ 2 ^^ i ' $ ^)] , (14) 

where IZij are the squared amplitudes, averaged over the incoming helicities and colors and multiplied 
by the flux factor l/(2s), for the NLO partonic subprocesses (gg — > 4>g,gq —> 4>q,qg — >• q(f>)'- 



3 Gji ag M 8 |^ gg (s,t,M)| s 
\/2 7T 2 s stw 



^(*l,«&rad) = £ V '^ 9V ;:' ^ , (15) 



„ ^ , G u a?M 4 s 2 + u 2 . . . , l2 

*»<•>••«> = - ^;„. (.+«)■« '^'■'■■ ) ' ■ (16) 

^ /5L ^ n G tf M 4 s 2 + i 2 . , l2 

^(*l,$rad) = ~ % ' \A qg {s,U,t)\\ (17) 

V 2 7T D S (S + tJ U 

where s = M 2 /(l - £), * = -(s/2) £ (1 + y) and u = -(s/2) £ (1 - y). 

The complete real matrix elements that enter the PDWHEG Sudakov form factor, eq. (jl]), read 

fl($l,*rad) = ^(^l,^rad)+^[i? OT (^l,^rad)+fi W (^l,^rad)] , (18) 

= , (19) 

where the functions R a b are the non- infrared-subtracted counterparts of eqs. ([8|)- (|10p . The probability 
for the emission of the first and hardest parton is described with the exact matrix element in all the 
phase space regions. 

Finally, the contribution of the qq —> 4>g channel is 

Rqq($l, $rad) = Kqq($l , $rad) Cqq, (20) 

with 



The functions A gg , A qg in eqs. (fT5j) - (fTT|) and A qg in eq. ([2Tj) depend on the particle content of the 
model considered, and will be defined in the following sections. 
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3 SM results 

The current public release of POWHEG [19] contains matrix elements evaluated in the HQET. It also 
gives the user the possibility of rescaling the term B(<&i) in eq. ([TJ by a normalization factor defined 
as the ratio between the exact Born contribution where the full dependence from the top and bottom 
masses is kept into account and the Born contribution evaluated in the HQET. In the following we 
describe the modifications we have introduced in the code to include the full fermion-mass dependence 
at the NLO and the effect of the two-loop EW corrections. 



3.1 Modifications in POWHEG 

The inclusion of the fermion-mass effects is achieved using for the functions H ie , T~L 2e , A gg , A qg , A qq 
the exact results instead of those computed in the HQET. For the Born term we have 



q=t,b 



where 



m 



2 



2 "(!-%,) \\r?{x q ) 



(22) 



Vq -M^ * 9 ~ ^1=^+1' (23) 

Tp = 1/2 is the matrix normalization factor of the fundamental representation of SU(N C ), and \ q is 
a normalization factor for the Higgs-quark coupling. In the SM case X q = 1 for both the top and the 
bottom quark. 

The form factor 1~L 21 contains the mass-dependent contribution of the two-loop virtual corrections, 
and can be cast in the following form: 



Tp ]T \ q (c F Gfjf R) (x q ) + C A Qf/f A \x q )\ + h.c. , (24) 

g=t,b V / 



n,2l 

rL 

g=t,b \ 

where Cp = (N 2 — 1)/(2N C ). Explicit analytic expressions for q^ Cr ^ Ca ^ given in terms of harmonic 

polylogarithms can be found in ref. [9]. It should be noticed that Q^ Cr ^ depends on the choice of 
renormalization scheme for the quark mass entering the one- loop part of the form factor. In ref. [9] 
expressions for with on-shell (OS) or MS parameters are presented. In our implementation 

we allow the choice among the OS, MS or DR renormalization schemes. 

Concerning the real emission contributions, we have for the gg —> Hg channel 

\A qq {s,t,u)\ 2 = \A 2 (s,t,u)\ 2 + \A 2 (u,s,t)\ 2 + \A 2 (t,u,s)\ 2 + \A 4 (s,t,u)\ 2 , (25) 

where the functions A 2 and A4 can be cast in the following form: 

A 2 (s,t,u) = T F Y \ q y 2 [b 1/2 (s q ,t q ,u q ) + b 1/2 (s q ,u q: t q )] , (26) 

q=t,b 

Ai(s,t,u) = Tp X q y 2 [ci/ 2 (s q ,t q ,u q ) + ci/ 2 (t q ,u q , s q ) + c 1/2 (u q , s q ,t q )] , (27) 

g=t,b 
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(28) 



Explicit expressions for the functions bi/ 2 {s q ,t q ,u q ) and c 1 / 2 (s q ,tq,u q ) are given in ref. [TO] . 
The function A qq (s, t, u) relevant for the qq — > Hg channel is given by 



A qq {s, t, u) = T F \ q y q d 1 / 2 {s q ,t q ,u q ) , 



(29) 



q=t,b 



and di/2(sq,t q ,Ug) can be found in ref. [TO]. Finally A qg (s,t, u) relevant for the qg — > Hg channel can 
be obtained from A qq -(s, t, u) via 

Agg(s, t, U) = Agg(t, S, u) . (30) 



The two-loop EW corrections are included as a factor (1 + 5 EW ) which multiplies the term B(<&i) 
in the first line of eq. (UJ. This choice follows from the current structure of PDWHEG where the qq —> Hg 
channel is kept apart, because it is not proportional to the Born cross section in the collinear limit. In 
the SM case, the values of the correction S EW as a function of the Higgs boson mass can be obtained 
from ref. [TH] . 

3.2 SM: numerical results 

In this section we present numerical results for the production of an on-shell Higgs boson in the SM. 
We focus our analysis on the inclusion of the exact quark-mass dependence in the NLO corrections 
and on the effect of the EW corrections. We also consider the effect of merging POWHEG with a PS. 
The results have been obtained for the LHC with center-of-mass energy of 7 TeV, using the following 
numerical values for the physical input parameters: = 1.16637 • 10~ 5 GeV -2 , mt = 172.5 GeV 
and m b = 4.75 GeV [3J. We have used the MSTW2008 [59] NLO set of PDF to describe to partonic 
content of the proton. In the code the value of a s (mz) is set accordingly to the choice made in the 
PDF set: in our case a s (mz) = 0.12018. When discussing the distributions in the Higgs transverse 
momentum pip, a cut pip 7 > 0.8 GeV has been enforced. The renormalization and the factorization 
scales have been set equal to the Higgs boson mass: hr = n, F = mn- 

In the left panel of figure UJ we plot the total Higgs production cross section, without acceptance 
cuts, in three different approximations: the dot-dashed line corresponds to the current public POWHEG 
implementation, in which the NLO-QCD corrections are computed in the HQET and are then rescaled 
with the exact Born cross section (which includes the full dependence on the top and bottom masses); 
the dashed line corresponds to the POWHEG implementation presented in this paper, where the complete 
NLO-QCD calculation is employed, i.e. the top and bottom contributions are treated exactly in the 
NLO corrections; the solid line also includes the effect of the EW corrections. 

In the right panel of figure [1] we plot the ratio between the Higgs production cross section obtained 
using our version of POWHEG and the one obtained using the current public version. The dashed line 
omits the effect of the EW corrections, while the solid line includes it. For ran < 160 GeV the 
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Figure 1 : Total cross section for SM Higgs production in gluon fusion at the LHC (7 TeV) , as a function of the 
Higgs mass, including different subsets of radiative corrections. We show the absolute predictions (left panel) 
and their ratio with respect to the current PDWHEG implementation (right panel). 



exact treatment of the quark masses results in an increase up to ~ 6% in the cross section, with a 
further increase (up to a combined ~ 10%) when the EW corrections are included. This effect is 
mainly due to the bottom-quark contribution, which is not negligible when the Higgs boson is light. 
For tuh > 180 GeV the quark-mass effects and the EW corrections have opposite sign, resulting in 
a ~ —2% correction in the Higgs mass range up to 2mt. Above the threshold for real top-quark 
production, where the approximation m< — > oo is not valid, the large corrections due to the quark- 
mass effects in the QCD contribution are partially screened by the EW corrections. In the rest of the 
section we discuss the kinematic distributions of a SM Higgs boson at NLO QCD. Since the effect of 
the EW corrections is very close to an overall rescaling of the total cross section, we neglect them in 
the following and focus on the effect of the QCD corrections. 

In figure[2]we show the rapidity distribution for a Higgs boson with mass run = 120 GeV. In the left 
panels, the dashed red lines correspond to the distributions obtained using the current implementation 
of P0WHEG, and the solid blue lines correspond to our implementation. In the right panels we plot the 
ratio between the distributions obtained with the two implementations. We compare two different 
approximations: in the upper panels we show the pure (i.e., fixed-order) NLO-QCD calculation, while 
in the lower panels we show the event distributions according to the basic formula of P0WHEG, eq. ([T]), 
which includes the effects of the Sudakov form factors, merged also with the PYTHIA QCD PS. As 
appears from the plots, the exact treatment of the quark masses results in a ~ 5% enhancement, 
uniformly distributed over the whole rapidity range. In the absence of acceptance cuts, the results 
obtained by combining P0WHEG and PYTHIA do not differ significantly from the pure NLO-QCD results. 
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Figure 2: Rapidity distribution for a SM Higgs with m# = 120 GeV. Left plots: in red (dashed) the current 
implementation of POWHEG, in which the NLO-QCD corrections are computed in the HQET and are rescaled 
by the LO cross section with exact top and bottom mass dependence; in blue (solid) the exact NLO-QCD 
corrections with full top and bottom mass dependence. The results are obtained at fixed-order NLO QCD 
(upper plots), or including the effects of the Sudakov form factor and of the PYTHIA PS (lower plots). Right 
plots: ratio of the result of exact calculation over the result of the current POWHEG implementation. 



In figure [3] we show the transverse momentum distribution for a Higgs boson of mass mn = 120 
GeV. In the left panels we compare the current POWHEG implementation with ours but, differently 
from figure [2j we show separately the pure NLO-QCD, POWHEG and POWHEG + PYTHIA calculations. 
In the right panels we plot the full results (blue, solid lines) and the ones obtained by introducing in 
POWHEG only the exact top-mass dependence (black, dashed lines), both normalized to the results of 
the current POWHEG implementation. 

The pure NLO-QCD calculation (upper panels) diverges for vanishing Higgs transverse momentum, 
although in figure [3] the divergence is masked by the lower cut on p!f . The plot on the right shows that 
for ptf < rat it is the inclusion of the bottom-quark contribution in the NLO corrections that gives 
rise to a positive correction with respect to the result obtained with the current POWHEG version, while 
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for pff > m t there is a substantial modification of the distribution, with a large negative correction 
driven by the use of the exact top-quark contribution. 

The inclusion of the effect of the Sudakov form factor (central panels) modifies the distribution 
in such a way that the latter vanishes in the limit pj. — > 0. As can be read from eq. ([1]), the 
probability of emitting the Higgs in association with a parton depends on the product A x R/B, 
where R is the squared matrix element for real-parton emission, B is the Born amplitude, and A is 
the Sudakov factor, which in turn is exponentially suppressed by R/B, see eq. In the current 
PDWHEG implementation, the emission probability is computed in terms of the ratio R(t,oo)/B(t,oo), 
where both R and B are evaluated in the HQETo For small ptf , the Sudakov factor with exact top and 
bottom mass dependence, A(i+6, exact), used in our implementation is smaller than the corresponding 
factor A(i, oo) used in the current POWHEG implementation, because R(t + b, exact) / B{t + b, exact) > 
R(t, oo)/B(t, oo). This inequality holds for two reasons: first, the p!f distribution is proportional to R, 
and R(t + b, exact) > R(t, oo) for pff < 200 GeV [36J; second, the inclusion of the bottom contribution 
reduces the LO cross section with respect to the result obtained in the HQET [6J. Thus, as shown 
in the right plot, for small pff the Sudakov factor suppresses the p^ distribution by almost 10% with 
respect to the result obtained in the current PDWHEG implementation. Since the emission probability 
is also directly proportional to the ratio R/B, starting from ptf ~ 30 GeV this factor prevails over the 
Sudakov factor, and the distribution with exact dependence on the quark masses becomes larger than 
the one in the current POWHEG implementation by up to ~ 15%. Finally, for p^ > mt the inclusion 
of the full top-mass dependence leads to a negative correction, similar to the one already observed in 
the pure NLO-QCD calculation. The inclusion of multiple gluon emission with the PYTHIA QCD-PS 
(lower panels) does not change dramatically - in the absence of acceptance cuts - the results obtained 
including only the hardest emission. 

Figure H] shows the transverse momentum distribution for a Higgs boson of mass mu = 500 GeV. 
The meaning of the different curves is the same as in figure [3j In the pure NLO-QCD calculation 
(upper panels), the effect of using a finite top mass in the NLO corrections is always negative and 
monotonically decreasing for increasing transverse momentum. Contrary to the light-Higgs case, the 
use of a finite top mass in the Sudakov form factor (central panels) yields an enhancement of the 
transverse momentum distribution at small values of pi^ . For p^ > 80 GeV the mass effects induce 
instead a dampening of the spectrum. As for the case of mu = 120 GeV, the introduction of the QCD 
PS (lower panels) does not modify substantially the distributions. As can be seen from the plots on 
the right, for such a large value of the Higgs mass the bottom contribution is highly suppressed. 

In order to illustrate the uncertainty arising from the use of different PS algorithms, in the left panel 
of figure [5] we compare the transverse momentum distributions obtained with our implementation of 
PDWHEG matched with either PYTHIA (red dashed line) or HERWIG (blue solid line), for a Higgs boson 
of mass mu = 120 GeV. In the right panel of figure Owe plot the ratio of the distribution obtained 
with P0WHEG+HERWIG over the one obtained with P0WHEG+PYTHIA. The HERWIG showering algorithm 

2 In the current POWHEG implementation only the total cross section is rescaled by the exact top+bottom LO result. 
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Figure 3: Transverse momentum distribution for a SM Higgs with inn = 120 GeV. Left plots: in red (dashed) 
the current P0WHEG implementation, in which the NLO-QCD corrections are computed in the HQET and are 
rescaled by the LO cross section with full top and bottom mass dependence; in blue (solid) the exact NLO-QCD 
corrections with full top and bottom mass dependence. The results are obtained at NLO QCD (upper plots), 
including the effects of the Sudakov form factor (middle plots), including also the effects of the PYTHIA QCD 
PS (lower plots). Right plots: the full NLO-QCD results (blue, solid) and the ones obtained by introducing 
in P0WHEG only the exact top-mass dependence (black, dashed), both normalized to the results of the current 
POWHEG implementation. 
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Figure 4: Transverse momentum distribution for a SM Higgs with inn = 500 GeV. Left plots: in red (dashed) 
the current P0WHEG implementation, in which the NLO-QCD corrections are computed in the HQET and are 
rescaled by the LO cross section with full top and bottom mass dependence; in blue (solid) the exact NLO-QCD 
corrections with full top and bottom mass dependence. The results are obtained at NLO QCD (upper plots), 
including the effects of the Sudakov form factor (middle plots), including also the effects of the PYTHIA QCD 
PS (lower plots). Right plots: the full NLO-QCD results (blue, solid) and the ones obtained by introducing 
in P0WHEG only the exact top-mass dependence (black, dashed), both normalized to the results of the current 
POWHEG implementation. 
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Figure 5: Transverse momentum distribution for a SM Higgs with ijih = 120 GeV. Left: comparison of the 
results obtained with P0WHEG+PYTHIA (red dashed line) and with P DWHEG+HERW I G (blue solid line). Right: ratio 
POWHEG+HERWIG over POWHEG+PYTHIA. 



yields a differential cross section larger by up to 10-20% with respect to PYTHIA in the region of small 
transverse momenta (p^ < 20 GeV), and smaller by 5-10% for large values of pif . We remark that this 
behavior is not specific to our implementation of P0WHEG. Indeed, the current implementation yields 
differences of similar size when matched with HERWIG instead of PYTHIA. 

In our SM analysis the top and bottom masses are renormalized in the on-shell scheme. Therefore, 
the dependence of the differential cross section on the renormalization and factorization scales does not 
differ significantly from the case of the HQET. As an example, for mu = 120 GeV and pif = 20 GeV, 
we find for our P0WHEG implementation an uncertainty band of [—16%, +21%] around the central value 
when the factorization and renormalization scales are varied in a range between 0.5 mu and 2 mu while 
keeping the ratio of the two scales in the range [0.5,2]. In the case of the HQET the corresponding 
uncertainty band is [—16%, +20%]. 

The phenomenological relevance of the finite quark mass effects can be established by comparing 
the size of the latter to the most accurate estimate of the theoretical uncertainty band for the Higgs 
transverse momentum distribution. This band can be found e.g. in ref. [33], where NNLO-QCD 
results in the HQET have been matched analytically with the effects of soft-gluon resummation. For 
mu = 165 GeV and 10 GeV < p^ < 60 GeV, the band width, for a variation of the factorization, 
renormalization and resummation scales in a range between 0.5 m// and 2m# (while keeping the ratio 
of any two scales in the range [0.5,2]), amounts to [—8%, +12%] with respect to the central value 
computed with all the scales set equal to mn- For the same value of mu and range of p!f , we find that 
the mass effects on the shape of the Higgs transverse momentum distribution range between —8% and 
+8%, i.e. they are of the same size as the uncertainty band of ref. |43j . 
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4 MSSM results 

We extended our implementation of the gluon-fusion Higgs production process in POWHEG to describe 
also the production of the neutral CP-even bosons of the MSSM, h and H. Differently from the case 
of the SM, the CP-even Higgs boson masses of the MSSM can be predicted in terms of the other 
parameters of the model. At the tree level it is customary to consider the pseudoscalar mass m A and 
tan/3 as free parameters which, together with m z , determine the masses m;, and mn of the CP-even 
Higgs bosons and their couplings to the quarks (as well as to the squarks). Radiative corrections, 
however, induce in the Higgs masses and couplings a dependence on all of the MSSM parameters (see, 
e.g., ref. [60]). When studying Higgs boson production in the MSSM, therefore, it is necessary to 
compute the entire spectrum of masses and couplings of the model in a consistent way, starting from 
a given set of input parameters. In our numerical analysis, we use the code Sof tSusy [61J to compute 
the MSSM spectrum starting from a set of running parameters expressed in the DR renormalization 
scheme. However, it is in principle possible to interface our calculation of the Higgs production cross 
section with other spectrum calculators (such as, e.g., FeynHiggs [62]) that adopt different choices of 
renormalization scheme for the input parameters. 



4.1 Modifications in POWHEG 

In order to describe the production of the neutral CP-even bosons of the MSSM, three modifications 
to the SM implementation of POWHEG are required: a rescaling of the Higgs-quark couplings in the top 
and bottom quark contributions, the introduction of all the contributions from diagrams involving 
superpartners, and, finally, a rearrangement of the EW corrections (which have been computed only 
in the SM case). 

For the production of the lightest scalar, h, the normalization factors for the Higgs-quark couplings 
entering the functions H u , H 21 -, A gg , A gg , A qq - in section [37X1 become 

cos a sin a 

X t = ~r-^ , *b = , (31) 

sin p cos p 

where a is the effective mixing angle that diagonalizes the radiatively corrected mass matrix in the CP- 
even Higgs sector. The corresponding normalization factors for the production of the heaviest scalar, 
H , can be obtained from eq. (|3ip through the replacements cos a — > sin a in At and — sin a — > cos a 
in X b . 

Diagrams with a squark running in the loop give an additional contribution to the one-loop form 
factor H 11 --. 

m u = 4T F J2^t y& [i + ife ln2 (*«■ )] > (32) 

"Iff. 

where the sum runs over the six flavors for the squarks qi and the two mass eigenstates (i = 1,2) for 
each flavor. The variables y q% and x q% are defined in analogy to y q and x q in eq. ([23]) . and the couplings 
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for the stop and sbottom mass eigenstates t\ and b\ to the lightest scalar h read 
sin a f 1 . _„ 1 
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sin j3 
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2 8 



1 + cos 29t I 1 sin 
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(33) 



1 + cos 20 b 



l--sin 



cos /3 



1 1 2 

• - sin 26b H n^b + tj m z sin 2/3 

2 8 



1 + cos 29b ( 1 — 77 sin 



(34) 



In the equations above \x is the higgsino mass parameter in the MSSM superpotential, A q (for q = t,b) 
are the soft SUSY-breaking Higgs-squark couplings, 9 q are the left-right squark mixing angles and 
9 W is the Weinberg angle. The couplings for the squark mass eigenstates t 2 and 62 can be obtained 
from the corresponding couplings for t\ and b\ through the replacements sin 29 q — > — sin 26 q and 
cos20q — > —cos29 q . The squark couplings to the heaviest scalar H can be obtained from the squark 
couplings to h through the replacements — sin a —> cos a and cos a —> sin a. 

The couplings for the up-type and down-type squarks of the first two generations can be obtained 
from the stop and sbottom couplings, respectively, by setting the quark mass and the squark mix- 
ing angle to zero. However, it can be seen from eqs. (|32|) - (|34|) that all contributions from the first 



two generations of squarks are suppressed by the ratio m|/m~. . Furthermore, there are significant 
cancellations among the contributions of the four squarks in each generation (indeed, the total con- 
tribution vanishes for degenerate squark masses). Therefore, in what follows we neglect the first two 
generations, and focus on the stop and sbottom contributions. 

Additional contributions to the two-loop form factor TL 2 ^ arise from diagrams with squarks and 
gluons, with four squarks, and with quarks, squarks and gluinos. In our POWHEG implementation we 
use the results of ref. [27] for the stop contributions, obtained in the limit of vanishing Higgs mass, 
and the results of ref. [3T] for the sbottom contributions, obtained via an asymptotic expansion in the 
superparticle masses. 

The functions A2 , A^ and A qq entering the real emission contributions in section 13.11 also receive 
additional contributions from diagrams with a squark running in the loop: 



AA 2 (s,t,u) 
AA 4 (s,t,u) 

AA q q(s,t,u) 



771 s 



t f^2^t y\ t c o (*« > ki > u m ) + c o ik 
t f^2~y- Vq< d o ( s q* » ki ' u ^ ) ' 



(35) 
(36) 
(37) 



where s qi ,t qi and u qi are defined in analogy to s q ,t q and u q in eq. 



Explicit expressions for the 



functions b (% , % ,u qi ), 



c o (% , % , Uqi ) and d (sq. , t qi , u qi ) are given in ref. [TO] . 
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Finally, we need to adapt the electroweak correction 5 EW to the case of the MSSM. Although a 
calculation of the contributions to 5 EW from diagrams involving superpartners is not currently available, 
we can obtain a partial estimate of the EW corrections in the MSSM by introducing in the SM result 
appropriate rescaling factors for the couplings of the Higgs boson. In particular, when the Higgs boson 
mass is below the threshold for real top production the EW correction in the SM is dominated by the 
contribution of two- loop diagrams involving light quarks, in which the Higgs boson couples to a gauge 
boson [T71 HU Q2] . We can therefore approximate the EW correction for the production of the lightest 
scalar h as 

(W « sin(/3 - a) ^2 [ReH u BeGg + Im^ w Im^ 2 /] , (38) 

where the one- loop form factor "H is computed in the MSSM (i.e., it contains both the quark and 
squark contributions) and the explicit expression for the two-loop EW light-fermion contribution 
can be found in ref. [19J. In the case of the production of the heaviest scalar H the factor sin(/3 — a), 
which rescales the Higgs-gauge boson couplings, must be replaced by cos(/3 — a). However, we recall 
that the approximation of including only the light-fermion contributions becomes less justified when 
m H > 2m t . 

4.2 MSSM: numerical results 

In this section we present numerical results for the production of the lightest CP-even Higgs boson, h, in 
a representative region of the MSSM parameter space. Events are generated with our implementation 
of POWHEG, then matched with the PYTHIA PS. We compute the total inclusive cross section, as well 
as the transverse momentum distribution, for the production of a light Higgs in gluon fusion, and we 
compare them with the corresponding quantities computed for a SM Higgs boson with the same mass. 
For the relevant soft SUSY-breaking parameters (and for ll) we choose 

m Q = mu = m D = 500 GeV , X t = 1250 GeV , M 3 = 2 M 2 = 4 Mi = 400 GeV, = 200 GeV, 

(39) 

where: rriQ, mu and mp are the soft SUSY-breaking mass terms for stop and sbottom squarks; 
Xt = At — \i cot f3 is the left-right mixing term in the stop mass matrix; Mi (for i = 1,2,3) are the soft 
SUSY-breaking gaugino masses. We consider the input parameters in eq. (|39[) as expressed in the DR 
renormalization scheme, at a reference scale Q of the order of the squark masses (in particular, we 
take Q = 500 GeV) . Our choice for Xt is modeled on the so-called "m™ 1 " scenario" , in which the stop- 
induced radiative corrections maximize the mass of the lightest scalar h, allowing it to satisfy the lower 
bounds from LEP even for relatively low values of the stop masses (for our choices of parameters the 
physical masses of the two stops are around 280 GeV and 660 GeV, respectively). We consider both 
signs for the parameter fj,, keeping in mind that, in our conventions, the tan /3-dependent corrections 
to the relation between the bottom mass and the bottom Yukawa coupling [63] enhance the Higgs 
couplings to bottom and sbottoms for fx < and suppress them for fi > 0. 
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We perform a scan on the parameters that determine the Higgs boson masses and mixing at tree 
level, m A and tan /?, varying them in the range Jf| 90 GeV < m A < 200 GeV and 2 < tan j3 < 50. For 
each value of tan /3 we derive the soft SUSY-breaking Higgs-stop coupling At from the condition on 
X t , then we fix the corresponding Higgs-sbottom coupling as Af, = A t . 

For each point in the parameter space, we use the code SoftSusy [61] to compute the physical 
(i.e., radiatively corrected) Higgs boson masses rrih and m#, and the effective Higgs mixing angle a. 
We obtain from SoftSusy @ also the MSSM running quark masses fht and fhf,, expressed in the DR 
scheme at the scale Q = 500 GeV. The running quark masses are used both in the calculation of the 
running stop and sbottom masses and mixing angles, and in the calculation of the top and bottom 
contributions to the form factors for Higgs boson production (the latter are computed using the DR 
results presented in refs. |27l I31j). As discussed in ref. [31] . the use of m&(Q) in the one-loop form 
factor for gluon fusion, induces potentially large contributions, enhanced either by tan/3 or by 
ln(mf/Q 2 ), in the two-loop form factor H 2e . We checked that our results are not significantly altered 
if we compute H ie in terms of the running bottom mass expressed at the lower scale Q = m^. 

In figure [6] we plot the ratio of the cross section for the production of the lightest scalar h in the 
MSSM over the cross section for the production of a SM Higgs boson with the same mass. For a 
consistent comparison, we adopt the DR scheme in both the MSSM and the SM calculations. The 
plot on the left is obtained with fi > 0, while the plot on the right is obtained with fx < 0. In order to 
interpret the plots, it is useful to recall that for small values of tua it is the heaviest scalar H that has 
SM-like couplings to fermions, while the coupling of h to top (bottom) quarks is suppressed (enhanced) 
by tan/3. In the lower-left region of the plots, with small tua and moderate tan/3, the enhancement 
of the bottom contribution does not compensate for the suppression of the top contribution, and 
the MSSM cross section is smaller than the corresponding SM cross section. On the other hand, 
for sufficiently large tan f3 (in the lower-right region of the plots) the enhancement of the bottom 
contribution prevails, and the MSSM cross section becomes larger than the corresponding SM cross 
section. For fi < the coupling of h to bottom quarks is further enhanced by the tan /3-dependent 
threshold corrections |63j . and the ratio between the MSSM and SM predictions can significantly 
exceed a factor of ten. Finally, for sufficiently large tua, i.e. when the couplings of h to quarks 
approach their SM values, the MSSM cross section is smaller than the SM cross section. 

It is interesting to note that for intermediate values of tua there is a band along which the two 
cross sections are similar to each other. Indeed, the observation of a scalar particle with cross section 
in agreement with the SM prediction does not necessarily imply that the Higgs boson is the SM one. 
However, as will be discussed below, a more detailed study of the Higgs kinematic distributions can 
help discriminate between the two models. 

To assess the genuine effect of the squark contributions (as opposed to the effect of the modifications 

3 We remark, however, that parts of the (rriA, tan/3) plane considered in our study have recently been excluded by 
searches for Higgs bosons decaying into tau pairs at the LHC albeit for different choices of the SUSY parameters. 
4 In the MSSM analysis we use M t = 173.1 GeV, m b (m b ) = 4.16 GeV and a s (m z ) = 0.1172 as inputs for SoftSusy. 
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Figure 6: Ratio of the total cross section for h production in the MSSM over the cross section for the production 
of a SM Higgs boson with the same mass. The plot on the left is for /i > while the plot on the right is for 
H < 0. 
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Figure 7: Ratio of the full cross section for h production in the MSSM over the approximated cross section 
computed with only quarks running in the loops. The plot on the left is for /i > while the plot on the right is 
for n < 0. 
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Figure 8: Left: ratio of the transverse momentum distribution for the lightest scalar h in the MSSM over the 
distribution for a SM Higgs boson with the same mass. Right: ratio of the corresponding shapes. 



in the Higgs-quark couplings), we plot in figure [7] the ratio of the full MSSM cross section for h 
production over the approximated MSSM cross section computed with only quarks running in the 
loops. As in figure O the plot on the left is obtained with \i > 0, while the plot on the right is obtained 
with n < 0. We observe that, in most of the considered region of the MSSM parameter space, the 
squark contributions reduce the total cross section. We identify three regions: i) for sufficiently large 
tan/3 and sufficiently small itia the squark contribution is modest, ranging between —10% and +5%; 
this region roughly coincides with the one in which the total MSSM cross section is dominated by the 
tan /3-enhanced bottom quark contribution, and is larger than the SM cross section; ii) a transition 
region, where the corrections rapidly become as large as —30%; this region coincides with the one in 
which the SM and MSSM cross sections are similar to each other; Hi) for sufficiently large uia the 
squark correction is almost constant, ranging between —40% and —30%; this region coincides with 
the one in which the MSSM cross section is smaller than the corresponding SM cross section. 

We now discuss the distribution of the transverse momentum p^ of a light scalar h, considering two 
distinct scenarios. First, we take a point in the MSSM parameter space (tha = 200 GeV, tan/3 = 10 
and n > 0) in which the coupling of h to the bottom quark is not particularly enhanced with respect 
to the SM value, so that the bottom contribution to the cross section is not particularly relevant. 
Because a light Higgs boson cannot resolve the top and squark vertices, unless we consider very large 
transverse momentum, we expect the form of the p\ distribution to be very similar to the one for 
a SM Higgs boson of equal mass, the two distributions just differing by a scaling factor related to 
the total cross section. This is illustrated in the left plot of figure [8j where we show the ratio of the 
transverse momentum distribution for h over the transverse momentum distribution for a SM Higgs 
boson of equal mass. In the right plot of figure [8] we show the ratio of the corresponding shapes, 
i.e. the distributions normalized to the corresponding cross sections. This ratio, as expected, is close 
to one in most of the pj, range. 
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Figure 9: Ratio of the transverse momentum distribution for the lightest scalar h in the MSSM over the 
distribution for a SM Higgs with the same mass. The six plots correspond to different choices of tua and tan /3 
for which the MSSM and SM predictions for the total cross section agree within 5%. The plots on the left are 
for \i > while the plots on the right are for fi < 0. 
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Figure 10: Ratio of the transverse momentum distribution for the lightest scalar h in the MSSM over the 
approximate distribution computed with only quarks running in the loops. The six plots correspond to different 
choices of tua and tan/3 for which the MSSM and SM predictions for the total cross section agree within 5%. 
The plots on the left are for \i > while the plots on the right are for fi < 0. 
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We then consider the opposite situation, namely when the coupling of h to the bottom quark is 
significantly enhanced. In this situation two tree-level channels, i.e. bb — > gh and bg —> bh, can also 
contribute to the production mechanism and influence the shape of the p T distribution |37j . Leaving 
a study of the effects of those additional channels to a future analysis, we will now illustrate how 
the kinematic distribution of the Higgs boson can help discriminate between the SM and the MSSM. 
The six plots in figure [9] correspond to different points in the (m,4,tan/3) plane characterized by the 
fact that the MSSM and SM predictions for the total cross section agree with each other within 5% 
(therefore, we are effectively comparing the shape of the transverse momentum distributions). The 
three plots on the left are obtained with [i > 0, while the three plots on the right are obtained with 
[i < 0. The Higgs boson masses corresponding to these points range between 114 and 122 GeV (i.e., 
a SM Higgs with the same mass as h would not yet be excluded by direct searches). Figure [9] shows 
that the region at small p T receives a positive correction with respect to the SM result for moderate 
values of tan/3. The correction decreases with increasing tan/3 and eventually becomes negative at 
large tan/3 for \x < 0. The region at large p T shows an opposite behavior with respect to tan/3. 

In figure [10] we show, for the same six points in the (m^,tan/3) plane as in figure the ratio 
of the p T distribution over the approximate distribution computed with only quarks running in the 
loops. Even though the light Higgs boson cannot resolve the squark loops, we see that the squark 
contributions can modify the form of the p T distribution, because of the interference with the bottom 
contribution. In particular, we observe that the squark contributions may yield a positive correction 
on the distribution at small pip, which turns negative for larger values of the transverse momentum. 
The negative correction becomes quite flat for p T > 100 GeV, and in the \i < case it reaches a —50% 
effect at very large p T . 



5 Conclusions 

We have presented a new implementation^] in the P0WHEG approach of the process of Higgs boson 
production via gluon fusion in the SM and in the MSSM. In the NLO-QCD contributions, we have 
retained the exact dependence on all the particle masses in the one-loop diagrams with real-parton 
emission and in the two- loop diagrams with quarks and gluons, whereas we have employed the ap- 
proximation of vanishing Higgs mass in the two-loop diagrams involving superpartners. We have also 
included the effects due to the two-loop EW corrections. 

The exact mass dependence at NLO QCD and its matching with the multiple gluon emission 
produced by the PYTHIA PS have important effects on the total and differential cross sections of the 
Higgs boson. In the SM, the exact dependence on the bottom-quark mass induces, for a light Higgs 
boson, a non-trivial distortion in the shape of the transverse momentum distribution in the small-p^ 
region. The effect is comparable in size with the current estimate of the theoretical uncertainty on 
this observable, which was derived in the HQET limit in ref. [43]. In the case of a heavy Higgs boson, 

5 A beta version of our P0WHEG code for gg — > <j> can be obtained upon request. 
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the role of the bottom quark is negligible, but the exact dependence on the top-quark mass yields, 
with respect to the HQET results, an increase of the transverse momentum distribution at small p!jZ , 
and a large negative correction at large p!f . 

In the MSSM, our code allows for a systematic study of the parameter space of the model in a 
realistic experimental setup. As an illustration, we considered representative choices in the MSSM 
parameter space, modeled on the so-called m™ ax scenario. We studied the role of the bottom diagrams 
and the impact of the inclusion of diagrams involving superpartners at NLO QCD, both on the 
total and on the differential cross sections. In the large-tan /3 regime, where the role of the bottom 
quark is very relevant, the differential distributions can receive large corrections, which cannot be 
described in the HQET approximation. A detailed study of the Higgs kinematic distributions could 
help discriminate between the SM and the MSSM, in case a scalar particle with a cross section 
compatible with the SM prediction is observed. 
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